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The hydro-dynamic expansion rate of quark gluon plasma (QGP) is evaluated and com- 
pared with the scattering rate of quarks and gluons within the system. Partonic scattering 
rates evaluated within the ambit of perturbative Quantum Choromodynamics (pQCD) are 
found to be smaller than the expansion rate evaluated with ideal equation of state (EoS) 
0^ • for the QGP. This indicate that during the space-time evolution the system remains out of 

equilibrium. Enhancement of pQCD cross sections and a more realistic EoS keep the partons 
closer to the equilibrium. 

I. INTRODUCTION 

The main aim of nuclear collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large 
Hadron Collider (LHC) energies is to create a thermalized system of quarks and gluons - mimicking 
the situation that prevailed after few microsecond of the big bang. The expansion time scale in 
the microsecond old universe is much larger than the interaction time scale of quarks and gluons. 
Consequently during the evolution the system maintains equilibrium. In the case of heavy ion 
collisions (HIC), however, these time scales have comparable magnitude. For HIC the two pertinent 
issues for the equilibration are (i) do the quarks and gluons achieve equilibrium and (ii) in case 
they do, can the equilibrium be maintained during the evolution ? or what kinds of EoS and 
interaction cross sections are required for the maintenance of equilibration. Some of these issues 
^-j- . will be addressed in the present work. 
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II. THERMAL EQUILIBRIUM IN HEAVY ION COLLISIONS 



To address the first issue the time evolution of average momentum {< pr >) of the quarks and 
gluons have been studied and compare with its equilibrium value. In relaxation time approximation 
it can be shown that the evolution of the < px > is governed by the equation: 

h ; 

. 5r . d<p T > _ <p T > - < p Tcq > 

dT Txelax 

Treiax is the relaxation time evaluated with pQCD cross section for a s = 0.3, PTeq is the value of 
Pt in equilibrium. The equilibrium time, r cq is obtained from the condition: < pt >— > PTeq as 
t - > T eq . It has been found that the value r cq depends weakly on initial value of < px >■ 

To address the second issue we assume that the quarks and gluons produced initially in thermal 
equilibrium and check whether it can maintain the equilibrium during the entire evolution processes 
by comparing their scattering rates with the expansion rate of the matter. The scattering time 
scales (T sca tt) of various partons are evaluated in pQCD and these time scales are compared with 
the expansion time scale (t cxp ). For maintenance of thermal equilibrium the following criteria 
should be satisfied: 

Tcxp ^ QT sca tt (2) 

where a ~ 0(1), is a constant. The criteria given in Eq. [5] is reverse to the one used to study 
the freeze-out of various species of particles during the evolution of the early universe [l[ ■ Similar 
condition is used in heavy ion collisions also to study freeze-out of hadrons [2| • Physically Eq. [2] 
means that the time taken by a partons to travel a distance of one mean free path by their thermal 



motion they have collectively receded from each other by less than one mean free path [2j]. The 
Tscatt is determined for parton types i by the expression 



^(JijVijUjj (3) 

where o\j is the total cross section for particles i and j, i?y is the relative velocity between i and j 
and rij is the density of the particle type j in the medium. 

The scattering time scales for partons have been calculated by taking into account the following 
leading order processes gg — ► gg, gg — » gg, g(g)g — ► g(g).g, qq — > 99, gg — > gg for light flavours and 
gluons Here g stands for light quarks and g denotes gluons. For evaluating r scatt for heavy 
quarks (Q) the pQCD processes are taken from The infrared divergence appearing in case of 
massless particle exchange in the i-channel has been shielded by Debye mass. 

The expansion time scale can be defined as: 



r" 1 - 1 &(T ' r) , 1 , 

cxp e(r,r) dr { ' 



where e(r, r) is the energy density at a proper time r and radial co-ordinate r. e(r, r) is obtained 
by solving the hydrodynamical equation: 

dpT^ = (5) 

with the assumption of boost invariance along longitudinal direction Q and cylindrical symmetry 
of the system In Eq. ([5|), T^ v = (e + P)u^u" — g^ v P is the energy momentum tensor, P is the 
pressure, u M denotes four velocity and g^ v stands for metric tensor. We consider a net baryon free 
QGP here, therefore the baryonic chemical potential (hb) is zero. 



III. INITIAL CONDITIONS AND EOS 



The expansion rates for RHIC and LHC energies have been calculated using the initial conditions, 
Ti = 400 MeV, n = 0.2 fm for RHIC which gives dN/dy - 1100 and Tj = 700 MeV, n = 0.08 fm 
for LHC giving dN/dy — 2100 The initial radial velocity has been taken as zero. Two sets of 
equation of state (EoS) have been used to study the sensitivity of the results on EoS. SET-I: In a 
first order phase transition scenario - we use the bag model EOS for the QGP phase and for the 
hadronic phase all the resonances with mass < 2.5 GeV have been considered 8]. and SET-II: The 
EOS is taken from lattice QCD (1QCD) calculations performed by the MILC collaboration Q. 



IV. RESULTS 



The time evolution of < px > obtained by solving Eq. [T] is depicted in Fig. [T] for a system 
of quarks and gluons. The results with pQCD interactions indicate that the quarks are unlikely 
thermalize both at RHIC and LHC energies. The possibility of gluon thermalization at LHC is 
very strong. The gluon thermalization time scale seems to comparable to the life time of the QGP 
at RHIC energies indicating that the gluons may remain out of equilibrium during the evolution. 
However, any non-perturbative effects may drive the system faster toward equilibrium. 

In Fig [5] the scattering time scale is contrasted with the expansion time scale for RHIC energy 
for two types of EoS mentioned above. Expansion rate for both the EoS is similar at early time 
but differ little at the late stage of the evolution. For the sake of comparison the expansion rate 
for the extreme case of free streaming is also displayed. The scattering rates are evaluated with 
pQCD cross sections (left panel). The condition for equilibration in Eq. [2] indicates that the gluons 
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FIG. 1: Evolution of partonic average momentum with time for a non-expanding system for RHIC and 
LHC energies. 



remains close to equilibrium, however the charm and bottom (not shown in the figure) quarks 
remain out of equilibrium during the entire evolution history. 

The analysis of the experimental data within the ambit of relativistic hydrodynamics suggest that 
the matter formed in Au+Au collisions at RHIC achieve thermalization. One possible reasons for 
the thermalization to occur is that the partons interact strongly after their formation. It is argued 
in [lfj | that the onset of thermalization in the system formed in heavy ion collisions at relativistic 
energies can not be achieved without non-perturbative effects. It has also been shown in [TTI ] that 
a large enhancement of the pQCD cross section is required for the reproduction of experimental 
data on elliptic flow at RHIC energies. Therefore, the pQCD cross sections used to derive the 
results shown in Fig. [5] (left panel) should include non-perturbative effects. To implement this we 
enhance the pQCD cross sections by a factor of 2. The resulting scattering time is compared with 
the expansion time in Fig. [5] (right panel). It is observed that the gluons are kept in equilibrium 
throughout the evolution, light quarks are closer to the equilibrium as compared to the heavy 
flavours. 

In Figs. [3] the results for LHC are displayed for the two time scales mentioned above for pQCD 
and enhanced cross sections. The expansion becomes faster at LHC than RHIC because of the 
higher internal pressure. As a consequence, it is interesting to note that the thermalization scenario 
at LHC does not differ drastically from RHIC. 



V. SUMMARY 



We have evaluated the rate of scattering of various partons in QGP and compared the scattering 
time scale with the expansion time scale to examine whether the evolving mater is in thermal 
equilibrium or not. It is found that gluons remain close to equilibrium throughout the evolution. 
Quarks gets equilibrated at the late stage of the evolution both for RHIC and LHC energies. 
Enhancement of pQCD cross sections and realistic 1QCD EoS drives the system toward equilibrium 
faster. 
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FIG. 2: Expansion versus scattering time scale for RHIC energy at r = 1 fm. 
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FIG. 3: Expansion versus scattering time scale for LHC energy at r = 1 fm. 
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